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Abstract 

We present new Monte Carlo results in non-compact lattice QED 
with staggered fermions down to too = 0.005. This extends our pre- 
vious investigations on the nature of the continuum limit of QED. 
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1 Introduction 



We report on the continuation of our investigation of lattice QED [1, 2]. The 
central question is whether the continuum limit of QED is trivial or not. We 
have been looking at non-compact QED with staggered fermions which has 
been found to have a second order chiral phase transition. 

In [1, 2] we have demonstrated that Monte Carlo data are consistent with 
the triviality picture. Two statements could be made. Direct evidence for 
triviality is that the renormalised charge seems to vanish when the lattice 
spacing is sent to zero. Indirect support for triviality is that the chiral con- 
densate a can be described by a mean field equation of state with logarithmic 
corrections. 

The author of [3] independently comes to the same conclusions when 
analysing the available Monte Carlo data. Both statements were also ar- 
rived at by employing Schwinger-Dyson methods both analytically [4] and 
numerically [5]. 

Nevertheless in [6, 7] it was argued that the chiral condensate should obey 
a power law equation of state with non mean field critical exponents. 

In order to try to resolve this controversy we have done new Monte Carlo 
runs on larger lattices down to a bare mass of m = 0.005. Additionally we 
directly measured the so called susceptibility ratio R [6]. It can be used to 
estimate the values of the critical coupling (3 C and the critical exponent 8 
without further assumptions. 

Using meson mass ratios as approximations of R the authors of [6] claim 
support for their ansatz of the equation of state. We have commented on 
the validity of these approximations in [2]. Here we will present direct mea- 
surements of R and relate these results to findings from the equations of 
state. 



2 Equations of State 

Our ansatz for the equation of state motivated from [1] is 



a 3 . a 



m ° = Al i — g/i / x +MP-Pc) , Pn , , (1) 

with parameters (3 C} p } q } Ai and A 2 . In our previous fit [1] we fixed q = 1 
and obtained (3 C = 0.187(1) and p = 0.61(2). The ansatz used in [6] reads 

m = A 5 a 8 + A b ((3 - /?> 6 (2) 

with parameters (3 C , 8 } A$ } Aj, and 6=1 (hyperscaling). The published result 
is 8 = 2.31 with (3 C = 0.205. 

With all parameters free both equations are numerically similar. From 
(1) one can derive effective values for 8 and 6 

8 eH = 3 + q/ log(l/<r), 6 eff = 1 + pi log(l/<r) . (3) 
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mean field EOS + log. corrections (1) 


(3 C q p 


X 2 /dof 


0.1887(3) 0.55(6) 0.53(2) 


5 


0.1869(1) = 1 0.65(1) 
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power law EOS (2) 


fi c 8 b 


X 2 /dof 


0.1885(3) 3.40(5) 1.28(1) 


8 


0.1954(2) 2.60(2) = 1 


30 


= 0.205 2.03(1) =1 


64 



Table 1: Fit results with MINUIT errors. 





chiral condensate at m 


= 0.01 


L 


f3 = 0.19 


fj = 0.20 


(3 = 0.21 
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0.1700(24) 


0.1097(18) 


0.0788(15) 


10 


0.1820(21) 


0.1244(18) 


0.0883(8) 


12 


0.1835(13) 


0.1322(10) 


0.0917(6) 


16 




0.1344(9) 


0.0965(5) 


oo 


0.1841(28) 


0.1348(13) 


0.0983(10) 



Table 2: Finite size dependence of a. The L = 10 data are taken from the 
second paper of ref. [6]. 

In Table 1 we have collected results from fitting (1) and (2) to Monte 
Carlo data including our new points at m = 0.005. We have done fits with 
all parameters free and also with fixed q = 1 and 6 = 1 respectively in order 
to compare with the results cited above. Values of \ 2 P er degree of freedom 
(x 2 /dof) are also given. The best fit is shown in Figure 1. 

Table 1 shows that our previous results [1] remain stable while the power 
law fit does not reproduce the parameters given in [6]. Using (2) with 6=1 
gives a relatively large x 2 /dof. Fixing (3 C to a value of 0.205 (the monopole 
percolation threshold) increases % 2 /dof considerably. 

Leaving all parameters free in (1) as well as in (2) leads to equivalent 
results. Positive values of q and p correspond to 8 e ^ > 3 and 6 e ff > 1 as seen 
from (3). These results are also compatible with our previous results [1] but 
contradict the claims in [7, 8]. 

Finite size effects might in principle influence the results. For the chiral 
condensate we have checked that these effects are small on our largest lattices. 
They are of the order of the statistical errors for the data we used in our fits 
(see examples in Table 2). 

To understand the finite volume dependence we have done calculations 
in the (3 — > oo limit. The leading dependence of the chiral condensate on the 
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Figure 1: Best lit of Table 1 to a data. 







linear lattice extension L was found to be [2] 

0"L - o"oo oc cr^exp(-3.23cr L iv) + 



(4) 



In Figure 2 we have plotted against the r.h.s. of (4). One obtains 
reasonably straight lines if ■ L is not too small. If we use (4) to extrapolate 
to L = oo we obtain the open points. Some extrapolated values are given in 
Table 2. We have also done fits to the extrapolated data. These give results 
that are consistent with those of Table 1. 



3 The Susceptibility Ratio 

The susceptibility ratio R is defined by 

d log cr da m 



R :-- 



d log m dm a 



(5) 
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Figure 2: Finite size scaling of a at m = 0.01. 



To independently measure R we employ the Ward identity da/dmo = C| ull (p = 
0) where C| ull (p = 0) is the full scalar correlation function at zero momentum: 

c s 11 = hll {(xx(x)xx(y)) - (xx(x))(xx(y))} ■ 

x,y 

The traces needed for the disconnected contributions were calculated with 
stochastic estimators. For each configuration we averaged over 20 vectors of 
random numbers. 

The plot of R versus m can be used to estimate (3 C and 8 [6]. In such 
a plot curves of constant (3 > (3 C bend upwards, curves of (3 < (3 C bend 
downwards and the curve for (3 = f3 c ends on the y-axis at 1/8. In Figure 3 
we show results of our measurements. From this plot one concludes that (3 C 
seems to lie near 0.19 and 8 seems to be close to the mean held value 3. 
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Figure 3: Measurements of R. 

4 Conclusions 

Our ansatz for the equation of state (1) describes the data well. If we also 
allow for hyperscaling violations in the power law ansatz (2) (i.e., keep b as a 
free parameter) we get consistent results with our ansatz. In particular the 
values of (3 C agree and 8 e ^ > 3. Demanding hyperscaling (b = 1) does not 
give such a good fit. We have checked that finite size corrections to these 
results are small. 

Our measurements of R are consistent with the results from the three 
reasonable fits to a. However, before we can draw definitive conclusions from 
the R data a closer look at the critical region is needed. The measurements 
are being extended to m = 0.005 and larger volumes. 

The main conclusion is that at present all Monte Carlo data are consistent 
with the hypothesis that QED has a trivial continuum limit. 
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